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NOTES 


INTRODUCTION 


This Physics Set of the Hewlett-Packard Computer Curriculum Series consists of a 
set of a Student Lab Book and a corresponding Teacher’s Advisor. It was designed 
to help meet the need for computer-oriented problems in physics, providing 
students an opportunity to use a computer as a problem solving tool within a par- 
ticular subject matter area. 


The materials are designed for flexible use as desired by the individual instructor. 
The material and exercises in this unit are intended as an “enrichment”’ experience 
in the field of geometrical optics. Fermat’s principle is used to discover the laws of 
reflection and refraction. Normally this requires differential calculus. However, 
with the speed of the computer, the same task requires only algebra and some 
trigonometry and therefore can be used at the introductory college level and in 
high school. The final topic covered is the method of ray tracing as applied to thin 
lenses. Since neither of the topics above is usually covered in the introductory 
texts, the use of this material will not compete with your text. Instead, it can be 
used to supplement and enrich in any fashion you choose. 


The degree of difficulty of the material is dependent upon the amount of assist- 
ance which you choose to provide. With no assistance, the better physics student 
should be challenged. However, given a good deal of assistance, any physics 
student should be able to work out the exercises with no great difficulty. The 
level of the material is determined by the assumption that students taking intro- 
ductory physics will be quite capable as a group. 


The Lab Book provides text material and programming exercises for the students, 
a sample program and advanced problems. The Teacher’s Advisor contains an ex- 
ample of a program to solve each exercise and a brief discussion of the important 
elements of the exercise. 


For best results, you should study all the solutions until you are certain you have a 
complete grasp of the general methods. This should be done before assigning the 
material to the class. Generally, the exercises are cumulative so that as techniques 
are developed they are used in subsequent exercises. Therefore, you will probably 
wish to proceed through the exercises in the order in which they are given. 


You will undoubtedly think of different programming methods or techniques as 
you study the example programs. Encourage the students to do the same. There 
are no approved solutions. All solutions are acceptable if they produce the correct 
results. At this level, there is no need for emphasis on the efficiency of a student’s 
program. The important question is, does it work? 
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REVIEW OF COORDINATE GEOMETRY 


Exercise 1 — Coordinates 


This is designed solely to provide practice in the location of points on a two- 
dimensional cartesian coordinate system. If your students do not need this prac- 
tice, go on to Exercise 2. The solution follows. 
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Exercise 2 — Distance Between Points 


The only question to be decided in this exercise is how to control the looping 
through the data. The program below uses a FOR loop. Other methods are pos- 
sible. In line 110, N is the number of line segments to be processed. Input of each 
segment takes place in line 130 with the computation in line 140. 


100 REM DISTANCE BETWEEN TwO POINTS 
110 READ N 

120 FOR I=1 TON 

130 READ X1sY1sX2sY2 

140 LET D=SQR((X2-X1)12+¢CYe-Y1)12) 
150 PRINT D 

160 NEXT I 

800 DATA 10 

801 DATA 1513536 

802 DATA -532s515-2 

803 DATA “3¢S5s-1le2sbe Bs heol 

S804 DATA 05905s031-4283 

805 DATA 03-6-732:030 

806 DATA *253-°35-3s-2 

807 DATA 09239-21492 —-27583--358 
808 DATA 051-41451-414s0 

809 DATA 2eBsle5s Hels 20] 

810 DATA 033+6s05s-2-4 

999 END 


RUN 


6240312 
7e2111 
8 «23286 
14283 
6732 
1¢41421 
169394 
19997 
7¢78267 
6 
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Exercise 3 — Rays Between Fixed Points 


The index of refraction is set equal to 1 in line 90. This is not required in this 
particular problem, but will be needed in subsequent ones. Hence, we establish 
the pattern here. In lines 100 through 130 the optical path P@ is computed for 
the first integer value of X, the value of X = 1 being assigned at the same time to 
the variable X@ which will be developed as the reflection point. A new optical 
path is computed in lines 150 through 170 each time through the loop. In lines 
180 through 190 the new optical path is compared to the smallest of the previous 
optical paths. If the new one is smaller, we update P@ and X@. 


90 LET Nei 

100 LET D1=SQR((1-1)124(00-4)12) 
110 LET D@=SQRC¢C11-1)9124+¢€6-0)12) 
120 LET PO=N*D14N*D2 

130 LET X0=1 

140 FOR X=2 TO 11 

150 LET D1=SQR¢CCX-1)124€0-4)72) 
160 LET D@=SQRCC11-X)124+(6-0)12) 
170 LET P1=N*D14N*«De2 

180 IF P1i>P0O THEN 210 

190 LET PO=P}) 

200 LET X0O=X 

210 NEAT X 

220 PRINT "REFLECTION AT X = "XO 
999 END 


RUN 


REFLECTION AT X = 5 
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Exercise 4 — Rays Between Changing Points 


This is a simple modification of the program in Exercise 3. The initial and final 
points are read into the program in line 120. The appropriate changes are then 
carried through the balance of the program. 


100 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
999 


RUN 


LET N=1 

INPUT X1s5Y1sX2sYe 

LET D1=¥Y1 

LET D2=SQR( (X2-X1)12+¥etred 
LET PO=N*D1+N*D2 

LET XO0=A41 

FOR X=X1it+1 TO Xe 

LET D1=SQR(CX-X1) 8e+Y1 12) 
LET D2@=SQR((X2-X) te+Yetead 
LET P1=N*D1+N*D2 

IF P1>PO THEN 240 

LET PO=P1 

LET XO=X 

NEXT X 

PRINT “REFLECTION AT X = “XO 
END 


702202310320 
REFLECTION AT X = 5 


Exercise 5 — Precision Computation 


S (in line 100) is the initial step size in the search. In line 110, M (number of 
significant digits in answer desired) is set. The search procedure in lines 120 through 
980 is the same as in Exercise 4. The only difference is found in line 180 which 
sets a counter C used to test against M. The value of the counter is tested in line 
290. If less than M the search limits are reset in lines 300 through 330. Then the 
program loops back to the search procedure. 


10 





100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
28 0 
290 
300 
310 
320 
330 
340 
350 
999 


RUN 


74 


LET S=1 

INPUT M 

LET N=1 

INPUT X1leYlsX2sY2 

LET D1=Y1 

LET D@=SQRC( CX2-X1)t2+yete2) 
LET PO=N* D1+N* D2 

LET X0=X1 

LET C=1 

LET A=X1+1 

LET B=X2 

FOR X=A TO B STEP §S 

LET D1=SQRCCX-X1)t2+Y1t2) 
LET De=SQRC(CXO-X) te4+yv2r2) 
LET P1=N*D1+N*« D2 

IF P1>P0O0 THEN 280 

LET PO=P1 

LET X0O=X 

NEXT X 

IF C >= M THEN 350 

LET A=X0-S 

LET B=X0+S 

LET S=S/10 

LET C=C+1 

GOTO 210 

PRINT “REFLECTION AT X = "XO 
END 


2702 103 1028 
REFLECTION AT X = 556101 


RUN 


24 


2-102 10910310 
REFLECTION AT X = -1+81608E-08 


RUN 


24 


759103109015 
REFLECTION AT X = 7+e00401 
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Exercise 6 — Angles of Incidence and Reflection 


The procedure in lines 100 through 340 is the same as in Exercise 5. Angles of 
incidence and reflection are computed in lines 350 and 360. These angles are con- 
verted to degrees and output in lines 370 and 380. 


100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
28 0 
290 
300 
310 
320 
330 
340 
350 
360 
370 
38 0 
999 


RUN 


24 


LET S=1 

INPUT M 

LET N=1 

INPUT X1sY1sXesYo 

LET D1l=Y1 

LET D@=SQRC CXe-XKX1) t2+yete2) 
LET PO=N*D1+4+Nx* D2 

LET X0O=X1 

LET C=1 

LET A=X1+1 

LET B=X2 

FOR X=A TO B STEP §S 

LET Dl=SQRC¢CX-X1) te4+y1t2) 
LET De=SQRC CX8-X) te4+yeare) 
LET P1=N*D1+N*D2 

IF P1>PO THEN 280 

LET PO=P1 

LET XO=X 

NEXT X 

IF C >= M THEN 350 

LET A=X0-S 

LET B=xX0+S 

LET S=S/10 

LET C=C+1 

GOTO 210 

LET I=ATNCCXO0-X1)/7Y1) 

LET R=ATNCCX2-X0) /Ye) 
PRINT “ANGLE OF INCIDENCE = '57.2957*I "DEGREES" 
PRINT “ANGLE OF REFLECTION ='57-2957*R"DEGREES" 
END 


2702 102203 5 
ANGLE OF INCIDENCE = 53-1335 DEGREES 
ANGLE OF REFLECTION = 53-1232 DEGREES 


12 
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Exercise 7 — Discovery 


The program in Exercise 6 is modified here to process several problems. The 
number of problems Z is input in line 95. The problem loop is established in lines 
96 and 400. Inside this loop, the details are the same as in Exercise 6. 


95 READ Z 

96 FOR J=1 TO Z 

100 LET S=1 

110 LET M=3 

120 LET N=1l 

130 READ XAlsYisxkesyYe 

140 LET D1=Y1 

150 LET De=SQRCCXe@-X1)t2e+vere) 
160 LET PO=N*D1+4N*D2 

170 LET X0=X1 

180 LET C=1 

190 LET A=X1+1 

200 LET B=Xe2 

210 FOR X=A TO B STEF S 
220 LET D1I=SQRCCXA1"-X)te+4+Y1172) 
230 LET De=SQRCCX2-X)128+Yer72) 
240 LET P1l=N*D1+4+N*D2 

250 IF P1I>PO THEN 280 

260 LET PO=P1 

270 LET XO=X 

280 NEAT A 

290 IF C >= M THEN 350 
300 LET A=X0-S 

310 LET B=xX0¢S 

320 LET S=S/10 

330 LET C=C+1 

340 GOTO 210 

350 LET ISATNCCX0-X1)/Y1) 
360 LET R=ATNCCX2-x0)/Y2) 
370 LET 12=57-2957*!I 

380 LET R=#57-2957*R 

390 PRINT X1isYilsxXesYesiIskR 
400 NEXT J 

800 DATA 5 

801 DATA -105953s854 

802 DATA 2510515210 

803 DATA 055210510 

804 DATA 20325230315 

805 DATA -203s63s-10310 

999 END 
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RUN 

-10 5 8 4 63-4577 63-4062 
2 10 15 10 33-0238 33-0238 
0 5 10 10 336636 33-7032 
20 25 30 15 14-0362 14-0363 
-20 6 -10 10 32-0051 3220055 


Fxercise 8 — Generalization 


In this exercise, the student is asked to “discover” the law of specular reflection. 
Note carefully that nothing in our development mentions the fact that the inci- 
dent and reflected rays lie in the same plane. This certainly has been implied since 
we have been working in a coordinate plane. As an optional exercise, students 
might be asked to repeat the development above in a three-dimensional coordinate 
system in which possible paths outside of a single plane could be investigated. Of 
course the student will discover that the rays do lie in a single plane. 


Exercise 9 — Corner Reflection 


This should be considered an optional exercise for the better student. Since two 
reflections are involved, a double loop is required to conduct the search. Trace the 
program carefully to see how the search is carried out. The object is to let the 
student discover the characteristics of a corner reflector. Have the students plot 
their results and discover that after two reflections, the ray is returning parallel to 
the initial ray. 


100 LET S=#i 

110 INPUT M 

120 LET N=! 

130 INPUT X1sYisxXesYe 
140 LET Dl=X1 


150 LET D2=¥Y1 
160 LET D3=SQR(X2t2+Yere) 
170 LET PO#N*(D1+D2+D3) 


180 LET X0=0 
190 LET YO=Y1 


200 LET C=1 


14 
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210 LET Al=Yie-l 

220 LET Bl=0 

230 LET A2=0 

240 LET Be=Xe2 

250 FOR Y=Al TO Bl STEP <5 
260 FOR X=Ae TO Be STEP S$ 
270 LET D1l=SQRCCY1-Y) tet+xXlt2) 
280 LET D@=SQRCYte+Xte) 
290 LET D3=SQRCY2t2+(Xe-X) te) 
300 LET Pl=N*(D1+4+De+D3) 
310 IF P1i>PO THEN 350 

320 LET PO=P1 

330 LET X0O=X 

340 LET YO=Y 

350 NEXT X 

360 NEXT Y 

370 IF C >= M THEN 450 

380 LET Al=YOtS 

390 LET Bl=YO-S 

400 LET A2&=X0-S 

410 LET Be=YOrS 

420 LET S=S/10 

430 LET C=C+} 

440 GOTO 250 

450 PRINT “REFLECTION AT" 
460 PRINT 0sYO3*AND"XO0s0 
999 END 


RUN 


23 
27531031035 


REFLECTION AT 
0 5e AND 5-« 0 


RUN 


23 
223953433 
REFLECTION AT 
0 22033 AND 1.75 0 
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Exercise 10 — Refraction With Precision Computation 


The modifications from the program in Exercise 5 follow from the discussion in 


the student’s Lab Book. 


100 LET S=1 

110 INPUT M 

120 INPUT N1sN2 

130 INPUT X1sY1,X2,Y2 

140 LET D1=Y1 

150 LET D@=SQR((X2-X1)124Y212) 
160 LET PO=N1*D1+N2*D2 

170 LET X0=X1 

180 LET C=1 

190 LET A=X1+1 

200 LET B=X2 

210 FOR X=A TO B STEP S$ 

220 LET Di=SQR((X-X1)12+Y112) 
230 LET D2=SQR((X2-X)12+Y212) 
240 LET P1l=N1*D1+N2*D2 

250 IF P1>PO THEN 280 

260 LET PO=P1 

270 LET X0=X 

280 NEXT X 

290 IF C >= M THEN 350 

300 LET A=X0-S 

310 LET B=x0+S 

320 LET S=S/10 

330 LET C=C+#1 

340 GOTO 210 

350 PRINT "REFRACTION AT X = "XO 
999 END 


16 





RUN 


74 

2151-33 

270210220s-10 

REFRACTION AT X = 12-693 


RUN 

74 

271¢3351 

70210s203-10 

REFRACTION AT X = 7-32201 


RUN 


24 

2151 

20310s205-10 
REFRACTION AT X = 10 


Exercise 11 — Refraction Between Different Media 


Following are runs of the program shown for Exercise 10 using the values given 
in Exercise 11. 


RUN 


24 
21246501433 
27-102102103;-10 


REFRACTION AT X = -1]-048 
RUN 

74 

71-3351-65 
2°1053103103s-10 
REFRACTION AT X = 2-085 
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RUN 

74 

7120051+-65 

27-102103s10s-10 

REFRACTION AT X = 428901 


RUN 

24 

21+6551-00 

2-102102105-10 

REFRACTION AT X = -4e27399 


RUN 


74 

2151 

27°-103103105-10 

REFRACTION AT X = -1-8&1608E-0& 


18 








Exercise 12 — Angles of Incidence and Refraction 


The changes in the program used for Exercise 10 are specified in the student’s Lab 


Book. 


100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 


290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
999 


LET S=1 

INPUT M 

INPUT Nl»sN2 

INPUT X1sY1sX@sY2 

LET D1l=Y1 

LET D2@=SQR(C(X2-X1)t2+Yere)d 
LET PO=N1*D14N2@*De2 

LET X0=X1 

LET C=] 

LET A=X1+41 

LET B=Xe 

FOR X=A TO B STEP S 

LET D1l=SQR((X-X1)?2+Y1 12) 

LET D@=SQR(CX2-X) te+Yere) 

LET P1l=N1*D14+N2@*D2 

IF P1>PO THEN 280 

LET FPO=P1 

LET XO=X 

NEXT X 

IF C >= M THEN 350 

LET A=X0<S 

LET B=XO0+S 

LET S=S/10 

LET C=C+1 

GOTO 210 

PRINT “REFRACTION AT X = "XO 
LET 12=57+¢29S5S7¥*ATNCCXO"-X1)/Y1) 
LET R=57+2957*ATNC CXE-X0)/Y2) 
PRINT “ANGLE OF INCIDENCE = “ABSCI1)"DEGREES” 
PRINT “ANGLE OF REFRACTION = “ABSCR)*DEGREES"” 
END 
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RUN 

24 

2124851233 
2-105105105-10 


REFRACTION AT X = -1-048 

ANGLE OF INCIDENCE = 41.8348 DEGREES 
ANGLE OF REFRACTION = 47-8504 DEGREES 
RUN 

24 


216335165 
27105103105 -10 


REFRACTION AT X = 2085 

ANGLE OF INCIDENCE = 50.393) DEGREES 
ANGLE OF REFRACTION = 38-3616 DEGREES 
RUN 

24 


212005165 
27-1051053105-10 


REFRACTION AT X = 4628901 _ 
ANGLE OF INCIDENCE = 55-0141 DEGREES 

ANGLE OF REFRACTION = 29.7306 DEGREES 

RUN 

24 


2126521-00 
2-103105105s-10 


REFRACTION AT X = -4.27399 

ANGLE OF INCIDENCE = 29.7955 DEGREES 
ANGLE OF REFRACTION = 54-9858 DEGREES 
RUN 

74 


2151 

2-105105105s-10 

REFRACTION AT X = -1-81608E-08 

ANGLE OF INCIDENCE = 44-9999 DEGREES 
ANGLE OF REFRACTION = 44-9999 DEGREES 


20 


Exercise 13 — Discovery 


The program below is a simple modification of the program used for Exercise 12. 
The modifications are specified in the Lab Book. 


100 LET S=i1 

110 INPUT M 

120 INPUT Nl»sNe 

130 INPUT X1sYisXesYe 

140 LET D1=Y1 

150 LET De@=SQR(CX2-X1)124yero) 

160 LET PO=N1*D1+4+Ne@«pDe 

170 LET X0O=X1 

180 LET Ce=12 

190 LET A=X1+41 

200 LET B=Xe 

210 FOR X=A TO B STEP §S 

220 LET D1i=SQRC(X-X1)124yY1172) 

230 LET De=SQR((X2-X)t28+vere) 

240 LET P1=N1*D1+N2*D2 

250 IF P1>PO THEN 280 

260 LET PO=P1 

270 LET X0=X 

280 NEXT X 

290 IF C >= M THEN 350 

300 LET A=X0-S 

310 LET B=X0+S 

320 LET S=S/10 

330 LET C=C+l 

340 GOTO 210 

350 PRINT “REFRACTION AT X = "XO 

360 LET I=ABSCATNCCXO0-X1)/Y1)) 

370 LET R=ABSCATNC (X@-X0)/Y2)) 

380 PRINT “ANGLE OF INCIDENCE = "57-2957*I"DEGREES"” 
390 PRINT “ANGLE OF REFRACTION= "57.2957*R"DEGREES" 
400 PRINT NI*IsN1*TANCI)sN1*COSCIDsNI*¥SINCI) 
410 PRINT N@*RsN2*TANCR)sN2*COSCR) sNE*SINCR) 
999 END 
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RUN 

73 

21¢4851-33 
27-105102105-10 
REFRACTION AT X = =-1-05 


ANGLE OF INCIDENCE = 41-8285 

ANGLE OF REFRACTION= 47-8555 
1-08047 1.3246 
12611087 1-46965 

RUN 

23 


214335165 
2-103105105s-10 
REFRACTION AT X = 2-08 


ANGLE OF INCIDENCE = 50-3815 

ANGLE OF REFRACTION= 38-3792 
161695 1-60664 
1-10524 123068 

RUN 

23 


712005165 
221021023103 -10 


REFRACTION AT X = 428 


ANGLE OF INCIDENCE = 54-9971 
ANGLE OF REFRACTION= 29-7695 
959882 1-428 
857302 ©9438 
RUN 
23 


21-6591-4-00 

2?°103103105-10 

REFRACTION AT X = -4-27999 

ANGLE OF INCIDENCE = 29-7695 

ANGLE OF REFRACTION= 54-9971 
©857302 °943801 
©959882 1-428 


22 


DEGREES 
DEGREES 
1-10261 


892431 


DEGREES 

DEGREES 
848104 
1629347 


DEGREES 

DEGREES 
°573616 
1.43225 


DEGREES 

DEGREES 
1.43225 
°573617 


987018 
°986137 


1.02451 
1.02443 


°©819124 
819246 


°819247 
819124 
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Exercise 14 — Generalization 


The student should come up with some version of Snell’s Law at this point. 


Exercise 15 — Double Refraction 


This program contains essentially the same search program as used in Exercise 10. 
The significant difference is that a double loop is required. 


100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
999 


LET S=1 

INPUT M 

INPUT Ni»sN2sN3 

INPUT X1sY1lsX2sY2 

INPUT ClsCe 

LET D1=Y1-Cl 

LET D2=C1i-Ce2 

LET D3=SQR( (X2-X1)124+CY2-C2)t2) 
LET PO=N1*Di1+N2*D2+N3*D3 

LET X8=X1 

LET X9=X1 

LET C=1 

LET Al=X1+1 

LET Bl=X2 

LET A2=X1 

LET Be=X2 

FOR U=Al1 TO B1 STEP S 

FOR V=A2 TO Be STEP S 

LET D1=SQRCCUHX1)1924¢CY1-C1)12) 
LET D2=SQR(CV-U) 824+(C1-C2)t2) 
LET D3=SQR((X2-V) 8t24CY1-Cldt2) 
LET PI=N1*D14+N2*D2e+N3*D3 

IF Pi>PO THEN 360 

LET PO=P1 

LET X8=U 

LET X9=V 

NEXT V 

NEXT U 

IF C >= M THEN 460 

LET Al=X8-S 

LET B1=X8+t+5 

LET A2=X9-S 

LET Be=x9t+S 

LET S=S/10 

LET C=C+l1 

GOTO 260 

PRINT “RAY PASSES THROUGH” 
PRINT X83C1”"AND *xX93C2 

END 
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RUN 
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2151-3351 

77103105105s-10 

2722-2 

RAY PASSES THROUGH 

-1.32 2 AND 1-33 -2 


Exercise 16 — A Complex Problem 


This program, like the one in Exercise 15, requires a double loop in the search. 
Essentially it is the same program, with a few geometrical differences. 


100 LET S=1 

110 INPUT M 

120 INPUT NisNe 

130 INPUT X1sY1sX2sY2 

140 LET Di=Y1 

150 LET D2=-X1 

160 LET D3=SQR(X2r1e+Yere2) 
170 LET PO=N1¥*(D1+D2)+N2*D3 
180 LET X0O=X1 

190 LET YO=O 

200 LET Cel 

210 LET Al=X1+1 

220 LET B1=0 

230 LET A2@=0 

240 LET Be=Ye2 

250 FOR X=Al TO Bl STEP §S 
260 FOR Y=A2@ TO Be STEP §S 
270 LET D1=SQR((X-X1)1284Y1192) 
280 LET De=SQR(Xr2+yYr2) 

290 LET D3=SQR(X2t24+C(ye-_Y) 12) 
300 LET P1=N1*CD14+D2)+4+N2*«D3 
310 IF P1>PO THEN 350 

320 LET PO=Pi 

330 LET X0=xX 

340 LET YO=Y 

350 NEXT Y 

360 NEXT X 

370 IF C >= M THEN 450 

380 LET Al=X0-S 

390 LET Bl=X0+S 

400 LET A2&=YO-S 
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410 LET Be=YOrS 

420 LET S=S/10 

430 LET C=C+1 

440 GOTO 250 

450 PRINT “RAY PASSES THROUGH 
460 PRINT X0,0"AND *O; YO 

999 END 

RUN 
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Fle1-5 


2210310510310 
RAY PASSES THROUGH 


“2-66 


0 AND 0 3°63 
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Ray Tracing 


RAY TRACING 


Exercise 17 — A Single Lens 
For the space: 
64 = 0, = —0.1 


y2 = y; +86; = 1 3 (20)(-0.2) = -3 


This ray is then carried through the lens: 
65 = 6, —y, /f = (-0.1) - (—3)/10 = +0.2 


y2 =yi = -3 


The output ray has @ = 0.2 and y = -3. 


Exercise 18 — Two Lenses 


The results of Exercise 17 are input into the space equations: 
> = 6, = 0.2 


y2 = y, + $6, = (-3) + (10)(0.2) = -1 


This ray is input into the lens equations to get the final result. 


6) = 6, —y,/f = (0.2) — (-1)/(-20) = -0.15 


y2 =y, =-!l 


The output ray has 6 = 0.15 and y = -l. 
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Exercise 19 — A Program for N Lenses 


Line 100 contains the description of the initial ray which enters the optical 
system. T stands for “Theta”. N] is the number of optical elements (either a 
space or a lens) in the system and is used to control the limit on the FOR loop 
opened in line 120. Inside the loop, the data for each element is read in, tested to 
see which equation is appropriate, and the ray then updated. After all the elements 
of the system have been processed the description of the output ray is typed out. 


100 READ YosT 

110 READ Nl 

120 FOR I=! TO NI 
130 READ N2sA 

140 IF N@=2 THEN 170 
150 LET Y=Y+A*T 
160 GOTO 180 

170 LET T=T-Y/A 
180 NEXT I 

200 PRINT YosT 

210 DATA Ose1 

220 DATA 4 

230 DATA 1,310 

240 DATA 2520 

250 DATA 1320 

260 DATA 2310 

999 END 


RUN 
2 -e15 
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Exercise 20 — Intersection of Rays 


Using the program from Exercise 19 the following results are obtained. 


210 DATA 130 
RUN 


0 --e95 


210 DATA ls-e} 
RUN 
-2 el 


Examining the results we see that the two rays differ in position by 2 units, have 
a relative angle between each other of 0.15, and are converging as they move away 
from the last element in the optical system. Thus the rays will intersect 2/0.15 = 
13.33 units to the right of the last element in the system. 


Exercise 21 — Description of Image 


The program from Exercise 19 is modified with new data and run with the results 
below: 


220 DATA @ 
230 DATA 1520 
240 DATA 2210 


210 DATA 150 


RUN 
i -el 


210 DATA 1s-e1 
RUN 


= 0 


The two output rays are two units apart, have a relative angle between them of 0.1 
and are converging. Therefore the image is real. The rays intersect 2/.1 = 20 units 
to the right | unit below the optical axis. Therefore the image is inverted, and the 
same size as the object. 
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OPTIONAL MATERIAL 


Exercise 22 — A Matrix Program 


The matrix X is used to hold the input ray. Each of the matrices corresponding to 
an optical element is read into the system and multiplied into M, the system 
matrix. Finally the output ray Y is the product of the system matrix M and the 
input ray X. 


100 DIM AC@s2I5sMC252),BC2s2)5XCEselI)5YC251) 
110 MAT READ X 

120 READ N 

130 MAT M=IDN 

140 FOR I[=1 TON 

150 MAT READ A 

160 MAT B=M*A 

170 MAT M=B 

180 NEXT I 

190 MAT Y=M*X 

200 MAT PRINT Y 

210 DATA e130 

220 DATA 4 

230 DATA 1Ls-els0al 
240 DATA 1520320,51! 

250 DATA 12--05033021 
260 DATA 15081051 

999 END 


RUN 
~e1l5 
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Exercise 23 — Description of Image 


This exercise is a very good test of understanding of geometry and all the material 
that has preceded. The rules are given below: 


Let the system matrix be described by four numbers as indicated. 


If a & c have the same signs, image is virtual. 

If a & c have opposite signs, image is real. 

Ifd > 0& b> O then the image is erect. 

Ifd>O& b< 0 & |d| > |bc/al| then the image is erect. 
Ifd>0&b< 0 & {d| < |bc/al then the image is inverted. 
If d<0& b < O then the image is inverted. 
Ifd<O0&b> 0 & {dl > |bc/al then the image is inverted. 
Ifd<O0&b> 0 & Idi < |bc/al then the image is erect. 

If |d — bc/al > 1 then the image is magnified. 

If Id — bc/al = 1 then the image is the same size. 


If |d — bc/al < 1 then the image is reduced in size. 
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